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Abstract 

Fluctuations in the Galactic synchrotron emission can be traced by the angular 
power spectrum of radio maps at low multipoles. At frequencies below few GHz, 
large-scale anisotropies are mainly induced by magnetic field turbulence, since non- 
thermal electrons radiating at these frequencies are uniformly distributed over the 
scales of magnetic field inhomogeneities. By performing an analysis of five radio 
maps, we extract constraints on turbulence spectral index and halo scale. Results 
favour a power spectrum significantly flatter than for 3D Kolmogorov-like turbu- 
lence, and a thin halo. This can be interpreted as an indication supporting non- 
conventional models of propagation of cosmic-ray particles in the Galaxy, or as a 
suggestion of a spectral-index break in the observed magnetic turbulence power 
spectrum. 
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1 Introduction 



One of the most important processes in tlie propagation of high-energy charged particles in 
the interstellar medium is their interaction with magnetic fields. The scattering of cosmic- 
ray (CR) particles on random hydromagnetic waves leads to their effective confinement in 
the Galaxy. This process of diffusion has a resonant character and it is mainly driven by the 
energy density associated to the random component of the magnetic field at the resonant 
wave number of the scattering which is given by the gyro-radius of the particle [1]. The 
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turbulence power spectrum at such scales and the extent of the turbulent region are key 
quantities in the description of CR propagation. 

The aforementioned scatterings change the energy of CR, reshaping their injection energy 
spectra into the ones we observe at Earth. The ratio of stable secondary to primary CR 
(where the reference ratio is often -B/C, given that Boron is almost entirely secondary) is 
a very useful diagnostic for such reshaping and in turn for the spectral index of turbulence 
power spectrum (for a recent review on CR, see, e.g., [2]), which is very often assumed to 
be a power-law P{k) ~ k~°' of the wavenumber k. Unstable secondaries are, on the other 
hand, a potentially very good tracer of spatial boundaries of the turbulent halo. Indeed, the 
ratio between stable and decaying isotopes depends on the CR confinement time, which is 
tightly related to the halo height (and the diffusion coefficient). 

In this paper, we aim to derive these two key quantities (turbulence spectral index and halo 
height) from a different perspective. We analyse the angular power spectrum of five radio 
maps. Their main component is given by the Galactic diffuse synchrotron emission. As we 
will discuss in Sec. 3, Galactic electrons involved in synchrotron radiation at frequencies < 
GHz are uniformly distributed over the scales of magnetic field inhomogeneities, given their 
much larger diffusive scale. Therefore the spatial spectrum of synchrotron intensity reflects 
the spectrum of magnetic fluctuations. By restricting to maps at low radio frequencies, we 
are confident to trace it avoiding contamination from electrons fluctuations, and to be in 
a regime where synchrotron emission dominates over other diffuse contributions (e.g., dust 
or free-free emissions) [3]. 

For CR, we consider a picture of fast particles scattering in a medium with weak hydro- 
magnetic turbulence and we will make the simplifying assumption of isotropic turbulence. 
The quasi-linear approximation leads to a diffusion coefficient D oc E"', with 7 = 4 — a [1], 
so its energy dependence is given by the same spectral index of hydromagnetic turbulence, 
which is derived by computing the angular spectrum of synchrotron emission and assuming 
a single power-law for the turbulence power-spectrum at all scales. 

The angular correlation, however, reflects statistics of turbulence only for sufficiently small 
angles, above which it follows a universal law ~ l/£ (see, e.g., [4]). Indeed, when the 
observation angle 6 for an observer located within an homogeneous turbulence region be- 
comes large, only points along of the hnes-of-sight very close to the observer are correlated 
while farthest points become uncorrelated. So as the angle increases the angular correlation 
w gets less and less contributions along line of sight, and follows the law w oc 6^^ which 
leads to Ci oc (see, e.g., Refs. [4,16]). The angle of transition between the two regimes 
is related to the ratio between the outer scale of turbulence L and the size of the turbulent 
region dmax- Provided an estimate for L (assumed to follow from scale of energy injection 
from supernovae remnant), dmax can be derived and represents an estimate for the height 
of the confinement halo in the Galaxy. 

Although all data considered in this paper come from single-dish surveys, next generation 
of radio interferometers with wide field of view and a primary beam (at low frequencies) 
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overlapping with the 'large' angular scales considered here, such as LOFAR [5], ASKAP [6], 
and SKA [7], can provide promising insight into the topic. Previous works on angular power 
spectrum of radio maps and its link with Galactic synchrotron emission include [8-17]. Their 
main focus has been on estimates of Galactic foreground and discussing its possible removal 
in CMB studies. In this paper, we instead focus the analysis on the extraction of turbulence 
properties and on their links to CR models. We also extend the determination of C^'s to 
few more maps. 

In Sec. 2, we describe datasets used and summarize the computation of C^. In Sec. 3, 
assumptions and model of angular power spectrum are discussed. Results on turbulence 
properties and their link to CR studies are presented in Sec. 4. Sec. 5 concludes. 



2 Data 



Numerous radio datasets concerning Galactic observations are currently available. We fol- 
lowed three main criteria to choose the maps used in this work: frequency < GHz, good sky 
coverage, and, possibly, good angular resolution. They include surveys of total intensity at 
22, 45, 408, 820, and 1420 MHz, and of polarization at 1420 MHz. Details are summarized 
in Table 1 and maps are shown in Fig.l in the HEALPix [25] format. ^ In order to avoid 
spurious projection effects, they have been obtained after regridding original maps into a 
much finer grid before filling pixels in the HEALPix tessellation scheme (with final linear 
size of pixels close to the original resolution of the survey) . 

The coefficients of the two-point angular correlation function are defined through the 
canonical relations T{d, (p) = Y^im (^iwXim{(^-, 0) and Q = l/(2£ -|- 1) Y}_i \aim^. 

The 'real' power spectrum coefficients can be estimated from measured C^"*' through 
Ci = (C^"^ — C""***^)/!^, which take into account effects of instrumental noise in (7""**^ = 
4:TT anoise/Npix (with anoise and Npix being, respectively, the rms temperature noise and 
number of pixels of the map) and of window function W — exp[—i{i + 1) (Jbeaml (with 
o'beam — FWHM/y/S In 2). The statistical uncertainty in the power spectrum can be esti- 
mated as 5Ci = {Ce + C^"''^/W) /^{2e + l)/2 Ai Uy, where Ai is the number of binned 
multipole and fsky is the fraction of sky observed. 

The C™"^ for maps listed in Table 1 are shown in central and right panels of Fig. 1. They 
have been computed using the Anafast tool of the HEALPix package [25]. We consider 
angular power spectra from all observed sky and from high-latitudes only [b > 20°, namely, 
cutting off emission from the Galactic plane). We then binned the interval £ — [5, 100] in 
20 bins of constant width in logarithmic scale. 



^ We acknowledge MPIfR's Survey Sampler [26] from which some of the datasets were down- 
loaded. 
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Frequency 


Beamwidth 


rms Noise 


Fraction of Sky 


p 

^max 


Survey 


MHz 


FWHM (°) 


K 


observed 






22 


1.1 X 1.7 


5000 


73% 


55 (30) 


DRAG [18] 


45 


3.6 


3500 


86% 


17 (n.c.) 


Guzman et al. [19] 


408 


0.85 


0.8 


100% 


100 (100) 


Haslam et al. [20] 


820 


0.97 


1.4 


51% 


55 f30) 


Dwingeloo [21] 


1420 (I) 


0.6 


0.02 


60% 


100 (100) 


Stockert [22] 


1420 (PI) 


0.6 


0.015 


66% 


100 (100) 


DRAG [23] 


1420 (PI) 


0.6 


0.012 


44% 


100 (100) 


Villa Elisa [24] 



Table 1 

Main parameters of surveys analysed in this work, imax = min(100, £n) is reported for all observed 
sky (and |6| > 20), where ^at is the multipole at which noise starts dominating, see text for details. 
At 1420 MHz, both total (I) and polarized (PI) intensity were considered. 



3 Method 

In the monochromatic approximation of synchrotron radiation, the energy of an electron 
emitting at frequency u is given hj E ^^\J ^ghz/ B^g GeV (where Vghz is the frequency 
in GHz and B^g is the magnetic field in ^G). It follows that, assuming a Galactic magnetic 
field of few microGauss [27], emissions at < GHz are mostly generated by electrons with 
energy significantly smaller than 10 GeV. The mean distance diffused by such electrons 
before losing most of their energy can be estimated as di = \/4:Dt [1], where D is the 
so-called diffusion coefficient and r is the the time-scale for the energy loss associated to 
radiative processes. In the energy range of interest, typical numbers for the Galaxy are 
D ~ Do - 3 • lO^^cm^s-i and r ~ 10^%E/GeV)-^ s [2], which translate into c^l > 4 
kpc for electrons with energy below 10 GeV. This order of magnitude estimate tells us 
that the scale of isotropy of the Galactic non-thermal electron population involved in the 
diffuse synchrotron emission detected in the radio maps is very large. Note that this is also 
because the storage time is far larger than CR electron injection rate which is given by 
supernovae explosions occurring as, roughly, few per century, which means that relative 
anisotropy generated by new energy injection is small compared to the total population of 
electrons (at low energy). For a numerical simulation of spatial distribution of Galactic CR 
electrons sec Sec. 4.1 in [38], which agrees with above arguments. CR electron fluctuations 
are negligible and does not contribute to the angular power spectrum (except for the very 
first multipoles < 5). 

Therefore fluctuations in the Galactic synchrotron emission can be totally ascribed to 
fluctuations in the magnetic fleld. Throughout the paper, we will make the assumption 
that the latter dominate the angular power spectrum for £ > 5 (6* < 30°). Fluctuations 
might be not all of stochastic nature, but with a contribution from the coherent magnetic 
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Fig. 1. ie/it Panel: Temperature maps (in K) of Table 1 in HEALPix format. At 1420 MHz only 
polarized intensity (obtained combining surveys of [23] and [24]) is shown. Horizontal lines indicate 
|6| = 20°. Middle Panel: l-cr regions of (binned) angular power spectrum Cj^""^ for maps in Table 
1. The normalization has been arbitrarily chosen (i.e., first bin is normalized to unity) to easy 
comparisons. The best-fit C^"^'' up to (-max is overlying (black solid curve). Note that it is obtained 
by fitting (C™°^ — C^°^'^^)/W with C^^^^ + C|, so it's not a direct fit to the shown observational 
C™°^, although strongly related. Right Panel: Same of middle panel, but cutting low latitudes 
(|6,„t| = 20°) data. 

field; on the other hand, in this case, fluctuations are expected to show up only at very 
large scales. Although synchrotron brightness does not directly trace the total magnetic 
field but only the component orthogonal to the line of sight, and a change in our viewing 
angle with respect to the smooth large-scale field (e.g., as our line of sight crosses a spiral 
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arm) can produce a sudden brightness variation, it is unlikely that such fluctuations can 
contribute to the angular correlation at small scales. The angular range discussed here can 
be considered to be descriptive of a truly turbulent regime. ^ 

We restrict to maps a,t u < GHz, disregarding, e.g., map of [28] at 2.3 GHz or WMAP [29] 
and Planck [30] maps, where the energy of emitting electrons would be larger and the 
associated mean diffusion length shorter, so fluctuations in the electron population could 
be sizable at smaller scales. This effect can also introduce a frequency dependence in the 
CiS, since the energy and, in turn, covered distances of emitting electrons are different at 
different frequency. On the contrary, for the frequency range considered here, we can flt all 
the datasets simultaneously, exactly because anisotropics are given by the same physical 
framework (i.e., magnetic turbulence) at any wavelength. We note that this implies that 
naive extrapolations of results found in the following Section to higher frequency (e.g., 
relevant for CMB studies) should be taken with a grain of salt, since they can be highly 
non-trivial. 

Note also that the non-thermal electrons should not be confused with the Galactic thermal 
electron population which instead shows turbulences (as we will mention in the next Sec- 
tion), but is not relevant for synchrotron emission; it can actually introduce anisotropics 
through free-free absorption of radio waves [31] but this effect is important only at < 10 
MHz. For the diffuse emission under investigation, synchrotron self- absorption is also not 
relevant. 

Both theory and numerical simulations support the idea of self-similarity of turbulence [1], 
which means that the spectrum of turbulence is typically predicted to follow a power-law 
A;^°, where a = 11/3 in the 3D Kolmogorov case. For this kind of spatial fluctuation 
spectrum, the angular power spectrum can be simply modeled [4,16] (for a > 1). Indeed 
turbulence scaling as /c^" implies Ci oc i^" in the small angle limit {9 < Oq) and Q oc 
in the large angle limt {9 > 9o), where 9 — n/i and the critical angle 9o can be 
estimated by ~ L/dmax with L being the outer scale of turbulence and dmax the boundary 
scale of the turbulent region. Although this assumes homogeneous turbulence and d„iax is, 
strictly speaking, the distance from the center to the farthest turbulence in a spherical 
scenario [16], one can take it as a rough estimate of the height of the turbulent diffusive 
halo; this is especially true when considering Q's from high-latitudes. Moreover, since, once 
again, inhomogeneity in synchrotron radiation are expected to be associated to magnetic 
inhomogeneity only, dmax should trace where magnetic turbulence density eventually drops, 
unless the turbulence has a peculiar inhomogeneity pattern. 

We model the angular power spectrum of turbulence as C^"''^ = C^pcf'^ /{Cf^ + C)"^) 
where Cf^ = Cii~^ and C^^"-* = C2^~" (namely, the two relevant scalings introduced above 
and in Sec. 1). Such form for C^"'"'' has been arbitrarily chosen to include scahngs at 



^ Throughout the paper, we make the assumption that small-scale magnetic field fluctuations 
trace the underlying interstellar turbulent motions, and so that measurements of the former 
(static quantity) provide information on the latter (dynamical phenomenon). 
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small ^ and £ " at large I. The exact shape at the transition between the two regimes 
affects only very mildly our results. 

At frequencies below 1 GHz, synchrotron contribution largely dominates the radio diffuse 
emission and other components as thermal bremsstrahlung and dust can be neglected [10]. 
We will however mainly refer to results obtained cutting off the Galactic plane such that 
this assumption can be more confidently trusted. 

At high multipoles the contribution of (extra-galactic) sources becomes important and needs 
to be included. On the other hand, in radio maps, it is typically found to be subdominant 
up to £ < 100 [13,15], so marginally relevant for our purposes. At £ ~ 100 (namely, in 
the multipole range where they can be possibly relevant), both observations [13,15] and 
theoretical models [32,33] suggest the contribution of sources to be approximately flat (the 
'Poisson noise' term dominates). Wc do not attempt to perform the delicate procedure of 
removing single sources, but rather we model the component, by adding an unknown flat 
contribution C| ~ = const at each frequency (and marginalizing over it when fitting the 
C,'s). 

Polarized intensity is potentially an even better tracer of synchrotron anisotropics than total 
intensity since it is not contaminated by thermal emissions. However, Faraday rotation 
effects around the source and along the line-of-sight (together with instrumental effects 
related to beamwidth and bandwidth of observations) can spoil the original power spectrum 
in a hardly predictable way. Faraday rotation is extremely important for radio data (while 
nearly irrelevant at higher frequencies) and multi- wavelength observations would be in order 
to single it out, but unfortunately they are not available. From Fig.l (bottom panel), one 
can see that C^'s from polarized intensity shows a slightly flatter trend than total intensity. 
If this is interpreted as purely due to Faraday depolarization effects, an estimate for scaling 
of the mean differential rotation of polarization angle A0 can be derived [14] . We obtain 
(sin A0/0)^ oc /Cj oc with n ~ 0.27 and 0.19 for all observed sky and high latitudes, 
respectively. Such a depolarization scaling (increasing with the angular scale) suggests an 
effect analogous to a beam depolarization (i.e., due to averaging of polarization vectors of 
different orientation); it is however difficult to separate different depolarization effects on 
C^'s, including that of, for example, fluctuations in the thermal electron population. Due 
to such uncertainties in extracting the intrinsic angular spectrum, we disregard polarized 
datasets in the fit of turbulence properties. 

Theoretical are thus given by C^"''* + C| which are described by four parameters 
(a, Ci, C2, Cg) as mentioned above; they have been determined by fitting to the expected 
(jjinap _ (jnmsey-^Y j-jj^^g^j ^j^g same way as C^*"^ in Fig. 1. In particular, we are in- 
terested in two parameters: a and Iq, which, in our model, do not depend on frequency. 
The transition scale io ~ ndmax/L between the two regimes and is defined by 
Cio"* = and can be easily expressed in terms of Ci and C2. The remaining two absolute 
normalizations of Cj"^^^ and C| can be seen as nuisance parameters (and can vary with 
frequency) . 
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Fig. 2. Contour plot of a and Iq for the mean likelihood and the 68% and 95% C.L. regions (black 
lines). The fit has been performed on all available data in all maps (left) and applying a cut for 
6 < 20 (right). See text for details. 

For each map, we consider multipoles up to Irnax = min(100, ^n) (reported in Table 1) 
where £jv is the multipole such that C™'*^ ~ C^°***^/IV. In this way, we conservatively avoid 
to aggressively extract information from multipoles where spurious effects and artifacts 
(related to projection, resolution, point sources, etc.) might have an impact on the results. ^ 
Indeed 're-projection' effect can affect only high multipoles, namely scales close to the 
angular resolution of the survey (unless data have been irreversibly binned to a significantly 
suboptimal resolution, which is not the case for the maps considered). 

Contour plots for a and Cq are shown in Fig. 2. They have been obtained considering a 
Gaussian likelihood function C = exp(— x^/2) and employing Bayesian statistics. ^ We 
assume flat priors for all parameters, so the posterior probability is basically given by the 
likelihood and the method corresponds to a maximum likelihood estimation where contours 
in Fig. 2 are given by Ax^ = ~ xlf^ where xlf is the of best-fit. This procedure does 
not imply parameters outside the confidence regions have necessarily a bad absolute x^? but 
such technique is more reliable than standard frequentist schemes when deriving preferred 
values of parameters for a given model, which is the aim of this work (for a recent review 
on the topic, see, e.g., [36]). 



^ The Haslam et al. map [20] is known to show scanning artifacts along lines of constant right 
ascension. We downloaded the map from NASA's LAMBDA website where original data were 
processed to mitigate baseline striping. In any case, Ci extracted from the original map or from 
a destriped version [34] shows negligible differences [15]. 

^ We acknowledge the use of sampling routines and an analysis tools of the CosmoMC package [35] 
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4 Results and Discussion 



By fitting Q's from all observed sky, we derive a = 3.1 ± 0.2 and = 10 ± 3 at 95% C.L.; 

considering high-latitudes data only, a = 2.91qj and io = 5ll\^ are obtained. The best-fit 
is found to be xlf — 63.6 (with 80 bins considered) and xlf — 65.3 (with 64 bins), 
respectively (we note also that fitting with a single power-law, one finds xlf ~ ^7.2 and 
xlf = 66.2; therefore such model, which however would lead to an even flatter spectrum, is 
disfavoured with respect to the one considered here) . The best fit of a is shghtly shallower 
at high-latitudes but compatible at 1-a level. One might expect io to be found slightly 
smaller at high latitudes which is indeed the case for best-fit values, although the two 
results are again compatible at 1-cr. We also check that best fits of fiat components C| 
are consistent with subdominant contributions in the whole multipole range i = [5, 100], 
which (together with the very mild variations of Ces from full-sky to Galactic plane cutoff) 
confirms assumptions about point sources discussed in previous Sections. 

We can confidently affirm that our results implies 7 = 4 — q;>0.7. As discussed above, the 
multipole of transition £q is related to the geometry of the turbulence (namely, corresponds 
to the angle above which turbulences start to be less and less correlated since lines of sight 
pass through turbulences of separated turbulent regions) and its fit provides the size of the 
turbulent halo in terms of the turbulence scale. We have dmax — " L/n < Q-^ ioopc ^P*" 
(where 100 pc is a typical scale for L as motivated below). 

The spectral index and halo scale of turbulence are extremely important parameters in 
the description of propagation of high-energy cosmic-ray particles in the Galaxy. Indeed, 
for isotropic turbulences (i.e., same assumption of above analysis) and in the quasi-linear 
approximation, the diffusion tensor simplifies to the scalar expression [1]: 

12nkresP{Kes) 3 (1 - 7) A^i"^ (^S^ ' 



where Vg — l/Kes — R/B is the gyroradius (with i?=particle rigidity), P{k) oc k'^"'^ is the 
turbulence power spectrum (again assumed to be a power- law), is a wavenumber where 
the random magnetic field assume the value SBl, and the power spectrum is normalized 
through /~P(fc) = 5i?i/47r. 

It is common practice to reduce Eq. 1 to D = /3Dq where /3 = v/c. The coefficient Dq is 
often assumed to be constant up to a vertical scale after which free escape of particles is 
considered. Or, similarly, Zh can indicate the scale at which the diffusion coefficient becomes 
large (e.g.. Do oc exp (z/zh) [37]). 

Recent derivations of spectral index 7 from cosmic-ray data can be split into two categories. 
If data are interpreted assuming a low convection wind (i.e., null or linearly increasing in the 
vertical direction), ratios of secondaries to primaries can be typically well- fitted by diffusion 
(plus reacceleration) models with Kolmogorov (7 = 1/3) [39] or Kraichnan (7 = 1/2) [40] 
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spectrum of turbulence. On the other hand, when convective effects are introduced, flatter 
turbulence spectra (7 = 0.75 — 0.85) are favoured [41-43]. More generally, CR models 
foreseeing processes other than diffusion to be responsible for the flattening oi B/C at 
high-energy (as e.g., convection or sources spectrum) allow large value of 7. 

The analysis in this paper supports this possibility. Note that our result is in agreement 
with previous works [8,9,12,13,15,16], all finding, with different uncertainties, a spectral 
index 7 ~ 1 for ^ < 100. 

This interpretation of our result is, however, based on a crucial and not obvious assumption. 
The relevant turbulence scale for scatterings of CR particle with hydromagnetic waves is the 
gyroradius of the particle (e.g., for GeV CR particles, Vg ~ 10^^ cm). This is much smaller 
than the scales of turbulence probed by our analysis. We have in mind a scenario where 
supernovae explosions and stellar winds inject turbulent energy at scales ~ L; then such 
energy is transferred through a direct cascade mechanism (with no dissipation) towards 
smaller scales. The resulting small-scale turbulence is expected to show the same spectral 
index of larger turbulence. However, one can argue that a straightforward extrapolation of 
the power spectrum down to the very small scales of CR scatterings can be a hazard. 

Observations of turbulent magnetic field at very small scales (i.e., of order of 10^^ cm) are 
lacking. On the contrary, fluctuations of thermal electron density in the nearby interstellar 
medium have been shown to be well approximated by a Kolmogorov spectrum [44], for 
scales ranging from 10^° cm to 10^° cm. Then it is commonly postulated that magnetic 
fields are frozen into the ionized interstellar medium and so magnetic turbulence should 
show the same power-spectrum [27]; however this possibility, although reasonable, needs to 
be proven. E.g., as a counter-example, at scales > 10^^ cm, observational evidences suggest 
a significantly fiatter spectrum for the magnetic component of the turbulence, as we will 
comment below. 

Turbulence observations most often rely on Faraday rotation measurements (RM). From 
RM only, on the other hand, it is not possible to separate fluctuations in electron density 
from magnetic field inhomogeneities (since RM — J B ■ Uedl). By combining emission mea- 
sure (EM = / nidi) with RM of polarized background radio sources, Ref. [45] estimated 
the magnetic turbulence at angular scales 1° < ^ < 10° (corresponding to scales between ~ 
10 pc and 100 pc and roughly overlapping with the scales £ < 100 considered in this work) 
finding a spectral shape significantly fiatter than in a Kolmogorov case, and, more precisely, 
a — 2.7 ± 0.1, which is compatible at the 2-a level with above results (and, in particular, 
with a is fitted from high- latitudes data, which is the most directly comparable case) . ^ At 
smaller scales, EM data were not available; using RM data only the turbulence was found 
to follow a similar power spectrum up to scales ~ 4 pc, below which a confident distinction 
between different scenarios is quite hard (but with Kolmogorov spectrum being perfectly 
viable down to the limit of the survey ~ 0.03 pc). Although the analysis in [45] rehes on 
some assumptions (e.g., same spectral index for electron and magnetic fluctuations) and 

^ The quantity derived in Ref. [45] was the structure function D (x. 6^. For small angles, the 
spectral index jjl is simply related to the spectral index of Q by )U = a — 2 [16]. 
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is limited to a small high-latitude region in the outer Galaxy (while significant variations 
of RM structure function across the sky have been reported [46-48]), the rough agreement 
with this work is very encouraging. ^ On larger scales, between 0.5 and 15 kpc (i.e., above 
the strictly turbulent domain which is < L, so possibly probing magnetic fluctuations of 
different origin) , Ref . [46] derived the magnetic spectrum from the ratio of RM and disper- 
sion measure (DM = / Ugdl) of pulsars in the Galactic disk, finding a nearly flat spectrum 
with q; - 2 ~ 0.4. 

Prom a theoretical point of view, it is not completely clear if a Kolmogorov-like spectrum 
should be expected. Many numerical simulations of 3D magneto-hydrodynamic turbulence 
have reported a shallower power spectrum, and the issue is currently under debate (for 
recent reviews see, e.g., [49,50] and references therein). 

To summarize this part, the state of the art about magnetic field turbulence is that from the 
observational data we have (including this work), which sample turbulence only at scales 
above few pc, the power spectrum is found to be much flatter than in the Kolmogorov case. 
At smaller scales, there are hints for a Kolmogorov power spectrum, but data which could 
lead to conclusive proofs or to conclusive exclusions of other models are lacking. 

So we can take two perspectives. We can assume results found in this work to hold all the 
way down to very small scales, and the above discussion on consequences for CR models 
follows. Conversely, we can assume a Kolmogorov power spectrum at small scales, and take 
the results as an evidence for a transition between two regimes, namely, a departure from 
self-similarity for the observed spectrum. For example, Ref. [45] considered the possibility 
that the Kolmogorov turbulence is in form of thin sheets which for large angular scales 
effectively appears as 2D turbulence. In this scenario, going from small to large scales, a 
transition from a 3D to a 2D turbulence regime would occur, where a 2D Kolmogorov power 
spectrum a = 8/3 is not far from our best- fit values. Note that this uncertainty (namely, 
a degeneracy between spectral and spatial properties of turbulence) is a general caveat in 
these analyses and is given by the fact that we are able to observe only 2D projections of 
underlying 3D spectra. 

For what concerns the halo height z^,, it is not easily constrained by GR data since they 
typically depends on the ratio Do/zh, rather than on only. The strongest bounds come 
from 'radioactive clocks', namely, unstable secondaries, as mentioned in the Introduction, 
with the most precise measurements being, at present, ^°Be/^Be. Neglecting local effects, 
they suggest a diffusive zone significantly thicker than the Galactic disc, i.e., Zh > 4 
kpc [39,41,42]; this is also supported by observations of significantly extended magnetic 
structures in external galaxies (for a recent analysis of a spiral disk galaxy, see, e.g., [51]. 

Another important observable for deriving magnetic halo properties is the Galactic total 

^ Note that rotation measures are sensitive to the line-of-sight component of the magnetic field in 
thermal gas, while synchrotron data probe the perpendicular component of the field illuminated 
by cosmic rays. However, under our assumption of isotropic and homogeneous turbulence, the two 
can be directly compared. 
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synchrotron intensity (see, e.g., [52]). However, in this case, the disentanglement of spa- 
tial properties of magnetic fields from cosmic-rays distribution properties, and of different 
magnetic components (regular and turbulent) themselves is very challenging. 

Although djnax cannot be immediately identified with Zh, it represents a rough estimate, 
in particular if extracted from observations at high-latitudes, see above. Such estimate 
crucially depends on assumptions about outer scale of turbulence. The dominant source 
of turbulence is believed to be supernova remnants (and superbubbles), injecting energy 
in the surrounding interstellar medium on scales L ~ 100 pc (for a review, see, e.g., [53]). 
Therefore, we infer a significantly thinner halo Zh 1 kpc, in apparent confiict with CR 
data. 

The disagreement could be alleviated if most SNRs are clustered with supershells of much 
larger extent than 100 pc (basically, ~ 4 kpc would require L ~ 1 kpc). However, 
although many HI supershell have been observed in the Galaxy [54] (and in external galaxies 
as well) extending up to kpc-scales (and presumably crated by SN explosions), it is unlikely 
that the mean characteristic outer scale of turbulence (i.e., the scale at which SNR energy 
is transferred to interstellar medium) is significantly larger than 100 pc, which follows from 
both simulation (e.g., [55]) and observational (e.g., discussion in [47]) arguments. 

On the other hand, the local environment can significantly affect abundances of radioactive 
species, since they travel short distances before decaying (of order of 100 pc). In particular, 
an underdensity of a similar size surrounding the Sun can exponentially increase the amount 
of decays (lowering the local ratio) [56] . Therefore the model of the local interstellar medium 
is crucial for the determination of Zh from unstable secondaries. Current observations [57] 
do not allow a very detailed description of the Local Bubble to precisely estimate this 
effect [2,42]. However, assuming a simple and reasonable underdensity of rh ~ 100 pc, 
bounds on Zh can be strongly relaxed and Zh '^1 kpc becomes fully viable [56,42]. Moreover, 
Ref. [41] found that the diffusive zone in models with large 7 is constrained to be rather 
small. 

A consistent (but not conventional) picture seems to emerge and deserves further investi- 
gation. 



5 Conclusion 

We have presented an analysis of the angular power spectrum in five radio maps at mul- 
tipoles 5 < £ < 100. At these scales, the CfS are a tracer of fiuctuations in the Galactic 
synchrotron emission, since contributions from sources or other diffuse components are sub- 
dominant. At frequency below few GHz, large-scale anisotropics in synchrotron radiation 
are mainly induced by magnetic turbulence, while the non-thermal electron population 
involved in the emission has a rather smooth spatial distribution. We assumed a sim- 
ple isotropic and power-law spectrum of turbulence P{k) ~ The theoretical angular 
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power spectrum thus scales as Ci oc before ~ ?>Q dmax J^V^' fi^i^d Ci oc 1'°" afterwards, 
where d^ax can be considered in first approximation as an estimate for the boundary of 
the turbulent diffusive halo. Results of fits are mildly dependent on latitudes. The spec- 
tral index is found to be a ~ 3, much flatter than Kolmogorov (a = 11/3) or Kraichnan 
(a = 7/2) cases often considered in the literature. Moreover, data favour a very thin halo, 
Zh<l kpc. 

If extrapolated to smaller scales, such Galactic turbulence scenario is consistent with 
cosmic-ray data, but pointing towards non-conventional models where convective winds 
play a significant role. On the contrary, assuming Kolmogorov turbulence to hold at scales 
comparable to the gyroradius of Galactic CR particles, our results can imply a non-trivial 
mechanism for transferring turbulent energy from large to small scales. 
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